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Abstract 

This work deals with the overdamped motion of a particle in a fluctuating one-dimensional periodic 
potential. If the potential has no inversion symmetry and its fluctuations are asymmetric and correlated 
in time, a net flow can be generated at finite temperatures. We present results for the stationary current 
for the case of a piecewise linear potential, especially for potentials being close to the case with inversion 
symmetry. The aim is to study the stationary current as a function of the potential. Depending on the 
form of the potential, the current changes sign once or even twice as a function of the correlation time of 
the potential fluctuations. To explain these current reversals, several mechanisms are proposed. Finally, 
we discuss to what extent the model is useful to understand the motion of biomolecular motors. 
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1 Introduction 



Protein motors play an important role in intracellular transport phenomena. In the cytoplasm of eucaryotic 
cells, certain macromolecules move along a complex network of periodic polymer tracks and transport 
organelles or vesicles containing chemicals. These ATP-hydrolysing macromolecules successively attach 
to and detach from the biopolymer while walking on it. Structural analysis of several motor proteins 
together with in vitro experiments have revealed that they undergo a cyclic sequence of conformational 
changes to convert the chemical energy from the ATP-hydrolysis into a unidirectional movement along the 
biopolymer. (see [p]] and the references therein). Although this mechanochemical conversion mechanism 
is not yet fully understood, there is no macroscopic thermal or chemical gradient in the cell medium 
that determines the direction of movement. Thus an interesting question arises: which properties of the 
system protein+track favour one direction of movement rather than the opposite? This is certainly a 
difficult matter: one knows motor proteins of the same family that move in different directions along the 
same track. 

Recently, some progress has been made in attempting to answer the above question. Apparently, 
the biasing of the catalytic cycle of the ATP hydrolysis is due to time-correlated chemical fluctuations 
induced by far-from-equilibrium concentrations of the reactants. To describe the motion of molecular 
motors, Magnasco considered the model of a Brownian particle moving in a one-dimensional ratchet-like 
potential driven by a symmetric stochastic force 0. A similar model has been used by Feynman et al. 

to illustrate the Second Law of Thermodynamics: Thermal fluctuations cannot produce a transport. 
(The example is originally due to Smoluchowski [Q].) But if the fluctuations of the force are time- 
correlated, the restrictions of the Second Law of Thermodynamics cease to apply and the particle shows 
a net movement that can even overcome the action of an external force: The asymmetric ratchet works as 
a mechanical rectifier. Both ingredients, the time correlation of the symmetric fluctuations and the lack 
of inversion symmetry of the potential yield a non vanishing current. 

Astumian and Bier [§] proposed a different picture for the motion of motor proteins. They studied the 
problem of a strongly damped particle moving in a ratchet-like potential that fluctuates dichotomously, 
i.e. between two different states [§]. The two different states of the potential model the electrostatic 
interaction of the motor protein with the periodic charge pattern of the biopolymer. The charge of the 
protein is changed when ATP binds to it. As a result, the potential changes as well. In this model, a net 
current is obtained at finite temperatures if the coloured potential fluctuations are asymmetric. Similar 
models have also been proposed by Prost et al. [§, 0). These authors argued that the two states correspond 
to the attached and the detached state of the motor protein. For a recent overview on this and related 
models we refer to [01, where the reader can find a large list of relevant literature on the subject. 

Although a simple dichotomous process contains some essential features of the movement of motor 
proteins, it can on no account provide a realistic description of the ATP hydrolysis. The main problem 
is that, for a ratchet-like potential, the direction of the current is fixed. It is therefore interesting to study 
other noise processes as well. Mielke recently developed a method that applies both to Magnasco's model 
of a fluctuating force ^ and to the model by Astumian and Bier QICQ. For a fluctuating sawtooth potential, 
he observed in several cases a current reversal for slightly different parameter sets of the multiplicative 
noise. Despite its simplicity, the model shows that small modifications of the motor protein suffice to 
make it move in the opposite direction. Other noise processes that can take not only two values as the 
dichotomous process can also be motivated from the biological situation of the motor protein. Typically, 
the motor protein undergoes several conformational changes, so that the interaction between the motor 
protein and the substrate must be described by several interaction potentials, not only two. 

In this paper, we apply the method developed by Mielke to treat the movement in a fluctuating poten- 
tial. Up to now, most efforts have concentrated on studying the influence of the noise parameters on the 
induced stationary current. Therefore calculations were performed taking the simplest asymmetric poten- 
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tial, a sawtooth ratchet. We rather focus on how the current varies when the geometry of the potential is 
changed. This question is of general interest, since one might expect new phenomena in a more complex 
potential. It is also of relevance for the question, whether or not these models can be used to construct a 
realistic picture of motor proteins. As we will show, the stationary current, which is the most important 
quantity to be calculated, depends essentially on the specific form of the potential. Therefore we conclude 
that a quantitative comparison between experimental findings and theoretical predictions from the simple 
model should be taken with care. A small change of the potential may change the picture quantitatively 
and even qualitatively. 

After defining the model, we derive a set of tridiagonal recursion relations to calculate the current and 
the stationary probability distribution. In Section 3 we take these relations as a starting point to discuss 
the limiting case of white noise and compute the current up to second order in the correlation time. In 
Section 4 we present exact numerical results for piecewise linear potentials. It turns out that direction 
and rate of the movement depend strongly on the details of the potential. Section 5 summarizes the main 
conclusions of the present work and discusses the relevance of our results for the motion of molecular 
motors. 



2 Definition of the model 

The overdamped Brownian motion of a particle in a fluctuating periodic potential can be described by a 
Langevin equation, 

dx , , dV (x) , — „ , , 

We have chosen the units so that the friction coefficient is unity. The second term on the right hand side 
describes the thermal noise, as usual one has 

(m=0 (£(05(O> = 8('-O- (2.2) 

The first term describes a fluctuating periodic potential. V(x) is periodic with period L. We restrict 
ourselves to the case where the amplitude z{t) of the potential V(x) fluctuates. This case has also been 
studied by Astumian and Bier [Q] and by Prost et al. [Qj. As mentioned above, they discussed only the 
case where z{t) takes two different values. We consider the more general case where z{t) is a Markov 
process described by a general Fokker-Planck equation 

d -^=M zP (z,t). (2.3) 

The operator M z has one eigenvalue Xo = 0. The corresponding right eigenf unction (j)o(z) is the stationary 
distribution of the Markov process z{t). The other eigenvalues are non-positive. The eigenvalue equation 
of M z is 

M z <\) n (z) = -Ktyn(z), (2.4) 

where (j)„(z) are the right eigenfunctions of M z . We assume that < X\ < X2 < Without loss of 

generality we take (z) > 0. The time-dependent autocorrelation function of z(t) satisfies 

(z(t) Z (t'))-(z) 2 oc e -^ (2.5) 

for large \t— The correlation time x is determined by the largest negative eigenvalue ofM z ,i.e. x = X^ 1 . 
We restrict ourselves to a class of processes for which 

Z$n(z) = ln,n+l$n+l(z)+Yn,n$n(z) + Yn,»-l<|>n-l (z), n>\. (2.6) 



3 



We assume that the stationary distribution (j)o(z) of z is normalized to unity. Due to the recursion re- 
lations ( ]2.6| ) the eigenfunctions (j)„(z) can be written as (j)„(z) = g n (z)$o(z) where g n (z) are orthogonal 
polynomials with respect to the weight function §o(z), 

J dzg„(z)gm(z)§o(z) = B n ,m- (2-7) 

This class of Markov processes is very general. It contains many processes that occur in typical situations 
such as the Ornstein-Uhlenbeck process, the dichotomous process, sums of dichotomous processes, and 
kangaroo processes. 



Having defined the stochastic process z(t) via the Fokker-Planck equation (2.3) it is natural to work 



with a Fokker-Planck equation for the problem described by ( |2.1| ). To do this one has to introduce a joint 
probability density p(x,z,t) for x = x(t) and z = z(t). The Fokker-Planck equation for p(x,z,t) can be 
written in the form 

fjp(V ''' ?) d f zf{x)-T^)p{x,z,t)+M z p(x,z,t). (2.8) 



dt dx \ dx 

In the following we will only consider static properties of our model. Therefore it is sufficient to cal 



culate the static joint probability density p(x,z). It satisfies (2.8) with vanishing left hand side. For the 
class of processes under consideration, it is suitable to expand p(x,z) in the complete system of right 
eigenfunctions § n (z) of M z , 

p(x,z)= Po (x)<\> (z) + Z(- l ) n Uz)p'n(x)- (2-9) 

n=l 



Using the equations ( |2.6[ ) one derives easily a set of recursion relations for p n (x). 

J = yo,of(x)po(x) - Tp' (x) -yo,if{x)p'dx), (2.10) 

lo,if(x)po(x) = X 1 p l (x)+y U if(x)p[(x) - Tp'l(x)-y h2 f{x)p' 2 (x), (2.11) 

yn-l,nf( x )Pn-l( x ) = KPn(x) +ln,nf(x)P n ( x ) ~ T Pn (•*) ~ 1nji+\f{x)p' n+ l (*) for » > L ( 2 -12) 

J is an integration constant which has a simple physical meaning, it is the static current. A slightly 
different derivation of these equations has been given in [|9|, [T^]. In general it is not possible to solve these 
equations, but they provide a good starting point for various approximations and solutions for special 
cases. One possibility is to solve (2.1C-2.12) using matrix continued fractions (see e.g. JTT|]). The 



continued fraction has to be truncated, it can then be evaluated numerically. Another possibility is to 
solve ( P-lPI - pH^I ) for small x perturbatively. This has been done to first order in [|lC|], and higher orders 
can be calculated straight forward. We give the results up to second order in the next section. 

In the special case where the potential V(x) is piecewise linear, the force f(x) is piecewise constant 
and ( )2 . 1 Q| — pTT^ ) become linear equations with constant coefficients that can be solved explicitly. The 
remaining task is to satisfy the continuity conditions for p n (x). For the simplest case, a sawtooth potential, 
this has been explained in detail in []|, 10]. A generalization to other piecewise linear potentials is straight 



forward. In Section [| we will present some results for a potential with three and four pieces. 

3 x-expansion 

As already discussed in JTo|], it is useful to construct a x-expansion for constant y n m in the case of a 
fluctuating potential. When studying a potential driven by a noisy force one often uses a x-expansion for 
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fixed D = To i A- In the case of a fluctuating potential, however, this would imply arbitrary large potential 
fluctuations, which is clearly unphysical. 

The x-expansion can be obtained using standard perturbation theory for linear operators [12]. But 
it is also possible to start directly from the recursion relations ( |2. 10| - |2~T2"| ). To obtain a x-expansion for 
constant y, hm , one uses 

Po(x) = p Q0 (x) + poi{x)t + Pq2 (x)t 2 + 0{t 3 ), (3.1) 
for po(x) and similarly 

/ = / + 7 lX + 7 2 x 2 + O(x 3 ). (3.2) 



J i = -jryo,o - yooV{x) - yooV(x) • (3-6) 



This ansatz implies that p n (x) = 0(l n ). Thus the lowest order terms can be obtained from ( 2.10 ), which 
yields 

Yo.qV(-v) 

poo(x)=Ce t (3.3) 

and Jo = 0. C is fixed using the normalization of poo(x). In the next order, the term containing p\ (x) in 
(2.10) becomes important. The equation can now be solved by means of a variation of the constant, 

p m {x)=d l \x)e t (3.4) 
For cW(;t) we obtain 

C(1)(X) = "7 Jo ( Jl+ yo,if(y){^f(y)Poo(y)\) e^dy + C™. (3.5) 

The constant can be fixed through Jjfpoi {x)dx = 0. For J\ one obtains using (0) = (L) 

fo,K,__ I Q L m 3 dx 

L Y0.0VW L 7 

J e t dx J e t dx 

The sign of J\ depends only on the sign of To.o/M- This means that to first order in x one never has a 
current reversal. Later in the discussion of our numerical results we will see that the current changes its 
sign as a function of x. Therefore it is interesting to see whether this behaviour can be obtained within 
the x-expansion. The calculations are straight forward, for additional details we refer to [10]. The final 
result for the second term J2 is 

h = CCy? u { (5to,o - Ti.i ) [ f(x) 2 f(x) dx + 7 °'° (Yl 7 ! 2 ~ Y °' 0) £ f(x) 5 dx) 

+C^ 1 {£d 1 Hx)f(x)f'(x)dx- £f(x)C^(x)dx}, (3.7) 

with C = [J L exp(yQ ! oV(x)/T) dx] 1 . Depending on Ti,i and on f(x) this expression may be positive or 
negative. This shows that for sufficiently large x one may have a current reversal. 

The x-expansion presented here has been constructed perturbatively. The n-th order is well defined 
if the 7i-th derivative of the potential V (x) exists. Thus the x-expansion is well defined only for analytic 
potentials, otherwise it breaks down. This happens if we take a piecewise linear potential. In that case 
already J2 is not defined. A similar observation is made in for a fluctuating force. In Jl3| ] it was shown 
how one can obtain an asymptotic x-expansion for that model. In general one obtains an expansion in 
y/x. But in contrast to the model with a fluctuating force, the first term J\ in the perturbative x-expansion 
calculated above is always correct in the present case. 
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4 Exact numerical results 



As already mentioned above, the recursion relations become differential equations with constant coeffi- 
cients if the force is piecewise constant. Let us divide the interval / = [0,L] into a set of disjoint intervals 
4 and let us assume that f(x) = fa if x G 4- An appropriate ansatz for p n {x) is 



Po(x) = Y, c r,kaol^i r) e a " >x + boMifx€l k , 

r 

Pn (x) = C r^n^ k ' * + b n,k if n > 1 , x € J k . 



(4.1) 
(4.2) 



Inserting these expression in the recursion relations one obtains a generalized eigenvalue problem for the 
coefficients and , 



/0 \ 

h 

^2 



V o 



a k 



J 



( tftu \ 
a\,k 



(4.3) 



(4.4) 



/ Yo,oAYo,i/* 
Yo,i A Yi,i A Yi,2/ t 

Yi,2/* '•• 

V 







(4.5) 



For details of the derivation and for an analytical solution of this problem for some special noise processes, 
we refer to 1 1C ]. In general, one can always truncate the matrices A& and at some large value of 



n and solve the eigenvalue problem numerically. In the following we will discuss results for sums of 
dichotomous processes. In this case it is possible to solve the eigenvalue problem analytically. The 
remaining task is to calculate the coefficients c r ^ and the current J. These quantities can be computed 
using the continuity of p n (x) for n > and p' n (x) f or n > 1 at the points where the force jumps from 
one value to another. These continuity conditions together with the normalization of p$ (x) yield a set of 
linear equations for the unknown coefficients c,-/t in the ansatz (4. 1 4/Q ) and the current J. The current can 
finally be expressed as a ratio of two determinants. A detailed description of this procedure has been given 
in JlO| ] for a sawtooth potential, i.e. for the case where f{x) takes two values. For a single dichotomous 
process the current has a fixed sign, but already for a sum of two or more dichotomous processes the 
current may change its sign as a function of the parameters, e.g. of X. The main motivation of the present 
work was to study the dependence of the current as a function of the potential. The simplest case is that 
of a force f(x) that takes three instead of two different values, i.e. 



fi tfx£h = [0,U) 
f{x) = { f 2 iixeh = [L l ,L 2 ) 
f 3 ifx e h = [L 2 ,L) 



(4.6) 



The parameters cannot be chosen independently. Using Lq=0, L3 =L and introducing the interval lengths 
Ak := Lk — Lfc-i, we have YX=\ fk^k = ^ due to the periodicity of V(x). When the potential is close to a 
simple sawtooth potential, for example if f\ and f% are negative, fa is positive and A3 is smaller than 
Ai + A2 one observes a behaviour that is similar to what has been obtained in [JT0[] . But if Ai + A2 is 
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close to A 3 new phenomena can be observed. If A\ + A 2 = A 3 and f\ = f 2 , the potential has inversion 
symmetry and consequently 7 = 0. In the following our goal will be to study what happens if one has a 
small deviation of the inversion symmetric case. For instance one can take a situation where Ai is much 
smaller than A2 and is much larger than |/2|. A typical form of the potential is shown in Figure 
1. We begin with the discussion of the behaviour of the system if z(t) is a simple dichotomous process 




Figure 1 : A slight symmetry broken three piece linear potential. 

that takes two values z\ and zi with equal probability. We will always assume that (z) = Yo,o > 0, i.e. 
Z1+Z2 > 0. Then the first term in the x-expansion yields a positive current. It turns out that, depending 
on the special choice of the parameters, one obtains a current reversal even in the case of a dichotomous 
noise process. An example is shown in Figure 2. Here we have zi =0 and Z2 > 0. Let us first assume 
that x is small. If z = zi, the particle is only subject to thermal fluctuations. When the potential switches 
to Z2, the probability to find the particle on the longest potential slope 3 is the same than to find it on 
any of the two other slopes. For not too large temperatures, the thermal perturbation in the equation of 
motion can be neglected. A simple calculation then shows that the time required to reach the minimum 
of the potential along two potential slopes is minimized if the slopes are equal. On average the particle 
moves to the right; as predicted by the x-expansion, a positive current flows. We now consider the 
opposite situation of sufficiently large x. Then one has a nearly adiabatic behaviour. But, as shown in 
[|To|], the current vanishes in the adiabatic limit. The first correction is of order 0(1 /x). Suppose that the 
temperature is not too large. If z = Z2, the particle will move towards the minimum of the potential and 
we can assume that it reaches a nearly static probability distribution oc exp(— zV(x)/T). This distribution 
is peaked near the minimum, but due to the asymmetry of the potential the probability to find the particle 
on the left hand side of the minimum will be larger than to find it on the right hand side. When z switches 
to zi = 0, the particle moves only due to the thermal noise. But since the initial probability distribution is 
asymmetric, the probability distribution at a finite time will be asymmetric as well. It is simply given by 
the convolution of the initial distribution with a Gaussian. Therefore, when z switches back to z = Z2, the 
probability to reach the next minimum on the left hand side will be larger than to reach the next minimum 
on the right hand side. This produces a net current to the left. Therefore the model shows a current 
reversal. This can be observed in Figure 2. For small values of the interval length Ai =L\, the region 
where the current is positive is very small, for larger values of L\ the region where the current is negative 
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Figure 2: The current as a function of log(x) for different potentials. The parameters of the noise are 
N = l,zi = and, zi = 2. The temperature is T = 0.04. Each potential has the common parameters 
L 2 = 0.5, L3 = 1, fa = —2, fa = 0.47 and different values of the interval length Ai = L\ and the slope fa. 
The corresponding values of fa are: L x = 0.01 -> fa w -0.439; L x = 0.03 -> fa « -0.372; Li = 0.05 -> 
/ 2 = -0.3; Li = 0.07 -► / 2 « -0.221. 



reduces rapidly. 

It is interesting to study what happens when one varies other parameters of the system, e.g. the 
temperature, or what happens for a more general noise process, which showed already a current reversal 
in the case of a sawtooth potential. In the following we present some typical results for a potential 
with the parameters fa = —2, fa = 0.47, L\ = 0.03, L 2 = 0.5, and L = 1. From these values one obtains 
fa rs —0.372 and AV = I/3A3I = 0.235. We have chosen these parameters since, as seen in Figure 2, 
they show the new current reversal very clearly, whereas for somewhat larger or smaller values of L\ the 
regions where the current is negative or positive becomes very small. 

Let us first discuss the current as a function of the temperature for a single dichotomous process. In 
Figure 3 we show results for the current as a function of log (x) for various temperatures. The dichotomous 
process takes the two values z\ = —0.3 and zi = 2.3. One observes that for a small temperature region 
around T = 0.027 the current is positive for small x, becomes negative when x increases and then becomes 
positive again. The third region is strongly enhanced when the temperature is smaller and it vanishes 
when the temperature becomes larger. To understand this new effect, we first consider the system at 
temperatures, for which the condition Zjmax(/(x)) 2 <, 2T <C z\vam{f{x)) 2 holds. Then one basically 
has the same situation as in the case where one of the coupling constants is zero and the other positive. 
When z(t) = zi, the motion of the particle in the potential is smeared out by the thermal noise. If we now 
slowly turn the temperature down, the particle will begin to feel the influence of the flat potential. Let 
us consider the nearly adiabatic case at sufficiently low temperatures. If z(t) = z 2 , the particle reaches a 
nearly stationary probability distribution. When the potential switches to z\, the distribution spreads out 
and the particle diffuses over several potential valleys. Since the potential slopes are now much flatter, 
the particle has not enough time to reach the stationary distribution before the potential fluctuates and 
we can by no means invoke the stationarity argument used above. Nevertheless, the dynamics is again 
determined by the geometry of the potential: Due to the steep slope 1, the particle surmounts the right 
slope 3 more easily than the other two ones (see Fig. 1). On average, one obtains a positive flow. 
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Figure 3: The current as a function of log(x) for a symmetric dichotomous process. The two values are 
Z\ = —0.3 and z 2 = 2.3. 



Instead of a single dichotomous process we now take a sum ofN—\ identical dichotomous processes, 
each again with the two values z\ and z%. The static distribution of z is given by 



Hz) 



i 

2^ 



k=Q 



§{z-(N-k)zi-kz2). 



and we have z\,2 = ±y, where y 
Y„, m are given by 



N 



{n+l)(N-n) 



Y 2 , Jn,n = (Z.) ■ 



(4.7) 



. After performing the rescaling j — > the coefficients 



(4.8) 



In [ [To| ] it was shown that for this class of processes the eigenvalue problem (43) can be solved analyti- 
cally. 

The question is now whether the effect observed for a single dichotomous process can be observed 
in this case as well, and if perhaps additional new effects occur. Next we present some results for a sum 
of two dichotomous processes. Let us first discuss the behaviour for large values of X. As in the case of 
a single dichotomous process we have a nearly adiabatic behaviour. But now, z takes three values 2zi, 
Z\ +Zi, and 2z2- For z\ < 0, z\ +Z2 > and small temperature the behaviour should be that of a single 
dichotomous process that takes the two values 2zi, and z.\ +Zi- Therefore, the sign of the current depends 
on the difference of z\ and zi, i.e. on j. If the condition \2z\ \ < z\ +Z2 holds, i.e. if y < \/2(z), it will be 
positive and otherwise, negative. Figures 4 and 5 confirm this qualitative discussion. 

At higher temperatures, the motion in the potential 2ziV(x) is smeared out and the system behaves 
again as in the case of a simple dichotomous process with zi = and Z2 > 0. From our discussion, we 
expect only one current reversal for all y > \[2 (z) . However, it turns out that for higher values of y and 
sufficiently large temperatures a weak positive current is observed when x becomes large (Figure 6). To 
understand the different behaviour in Figure. 5 and 6, first note that for y = 1 .7 the quantities \2zi \ = 1 .4 
and z\ +Z2 = 1 have a similar numerical value, whereas for y = 2.5 one has 2z\ ~ —3.4 and z\ +Z2 = 1- In 
the first case the motion for both values z{t) is determined mainly by the thermal noise when T rises; but 
in the second, only the motion in the potential that corresponds to the value z(t) = (z) will be smoothed 
out while the motion in the potentials for the other two values is not influenced so much by the thermal 
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Figure 4: The current as a function of log(x) for a sum of two dichotomous processes at different temper- 
atures. The average of the coupling constant and the noise strength are (z) = 1 and y= 1. 
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Figure 5: The current as a function of log(x) for a sum of two dichotomous processes at different temper- 
atures. The parameters of the noise are (z) = 1 and y= 1.7. 
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Figure 6: The current as a function of log(x) for a sum of two dichotomous processes at different temper- 
atures. Here we have set (z) = 1 and y = 2.5. 



noise. For large x and not too high temperatures only the values with small \z{t)\ are relevant, i.e. 2z\ and 
Z\ +Z2- Thus the system behaves as in the case of a single process where one of the coupling constants is 
zero and the other is negative. This yields a small positive current. 

Let us briefly discuss the behaviour of the current for a sum of N identical symmetric dichotomous 
processes. Figure 7 shows some plots of the current as a function of x for different values of N. For 




-0.005 



Figure 7: The current as a function of log(x) for a sum of dichotomous processes. N takes the values 
1,2,3,4 and 5. The other parameters of the noise read (z) = 1 and y = 1.3; the temperature is T = 0.027. 



small values of the correlation time (< 10 2 ), all curves merge into a single curve. This is a generic 



feature, which holds also for general potentials V(x). The reason is that, according to (3.6) and (3.7), the 
coefficients J\ and J% of the x-expansion are the same for all N, since the quantities 70,0, Yo,i an d Yi,i do 
not depend on N. In contrast, the behaviour of each curve is different near the adiabatic limit; for Af > 1 , 
the effect responsible for a positive current disappears. 
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The fact that all plots lie so close to each other suggests that they rapidly converge to a single curve 
for N — > °°. In this limit, the sum of dichotomous processes yields a well known Gaussian process, the so- 
called Ornstein-Uhlenbeck process [14]. In this case the potential, rather than jumping between different 
states, suffers a continuous distortion. One therefore expects that some mechanisms that lead to a current 
reversal cancel. Nevertheless, our results show that the current changes its sign once, in contrast to the 
case of a sawtooth potential, where the sign of the current is fixed [10]. This behaviour is also predicted 
by the x-expansion, which leads to a current reversal for potentials of the form discussed above and an 
Ornstein-Uhlenbeck process. 

The validity of the expansion around the white noise limit can be tested with the help of the exact 
numerical results. In Figure 8 two exact plots of the current for N = 1 and N = 2 and the first order term 
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0.0000 
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-0.0004 
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log(x) 

Figure 8: Two plots of the current as a function of log(x) computed for a simple dichotomous process 
(solid line) and a sum of two processes (dashed line) are shown together with the first order term (dotted 
line) of the x-expansion. Here we have set (z) = 1 and y = 0.6. The temperature is T = 0.04. 

of the x expansion are depicted. The quantitative agreement with the exact curves is remarkably good for 
values of x below 10~ 4 . 

To conclude this section let us discuss briefly what happens for more complicated potentials. It is 
straight forward to do the same calculations as above for e.g. a piecewise linear potential with four (or 
more) pieces. We have done several calculations for a potential with four pieces. The results for small 
and large values of X are similar to the results reported for the potential with three pieces. For small x a 
current reversal may occur and can be predicted by the x-expansion. For large x, in the nearly adiabatic 
limit, the current is positive or negative, depending on the temperature and the noise parameters. The 
mechanism we discussed in that case for the potential with three pieces can be carried over to the more 
general case. For intermediate values of x, the current shows a more complicated behaviour. Depending 
on T and y, one observes additional minima and maxima for J as a function of x. On the other hand, we 
never observed more than two current reversals for J as a function of x. 




5 Conclusions 

The main result of this paper is that a new kind of current reversal occurs for a Brownian particle in a 
fluctuating potential, if the potential is not a simple sawtooth potential. For a saw-tooth potential a current 
reversal has only been observed for more complicated noise processes, but not for a simple dichotomous 
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process or for the Ornstein-Uhlenbeck process [§, |(| |lC[ ]. In a more complicated potential a current 
reversal even occurs for the dichotomous process and the Ornstein-Uhlenbeck process. It can even be 
predicted by a simple x-expansion carried out to second order. For a general potential we have derived 
the first two non-vanishing terms of the expansion around the white noise limit and we have shown that, 
depending on the noise parameters and the shape of the potential, they may have different signs, thus 
leading to a current reversal. 

The stationary current depends strongly both on the statistics of the coloured fluctuations and the 
details of the potential. Depending on the range of x and T , the current has a different direction and 
magnitude. 

We have focused on the special case of a piecewise linear potential. Our numerical results show 
that small deviations of the inversion symmetric case change the qualitative behaviour of the system 
dramatically. If one takes for instance a slight asymmetric three piece linear potential, the current may 
change its sign as a function of x more than once, even in the case of a simple dichotomous process. The 
number of current reversals depends on the relative sign of the coupling constants z\ and z.i- If both are 
positive or one is positive and the other one is zero, only one current reversal is observed, whereas in the 
case z\ < and z\ +Z2 > a new low temperature effect takes place and provokes an additional current 
reversal for large x. In this situation, the induced current tends to a finite value when T — > 0. 

Depending on the time scale of the fluctuations, the asymmetry of the potential acts on the particle in 
a different way: In the nearly adiabatic limit, the thermally driven particle surmounts one of the adjacent 
potential barriers more easily than the other one, whereas for small x (and not too high temperatures), 
the different drift times down the three potential slopes favour one direction of movement. Both effects 
are purely dynamical, it is not possible to explain the behaviour using transition rates calculated in the 
adiabatic limit. 

For a sum of two dichotomous processes, the system may be described by an effective dichotomous 
process if x is sufficiently large. The behaviour depends on the sign of the average of the coupling 
constant for the effective process, i.e. on (z) and y. When T rises, some values of z(t) will be washed out 
by the thermal noise. Therefore, the choice of the coupling constants of the effective process depends on 
temperature. 

In the general case of a sum of N dichotomous processes, we have seen that the current does not 
depend on the noise details for sufficiently small x. Whereas in the case of a sawtooth potential the 
current changes sign only if N is odd [10], for a proper form of the three piece potential one observes at 
least one current reversal for all values of N. 

Since our numerical results have been obtained only for a piecewise linear potential, we would like 
to emphasize that our results hold for general potentials as well. As mentioned above the current reversal 
in the case of a dichotomous process or in the case of the Ornstein-Uhlenbeck process can be predicted 
by a second order x-expansion for any given potential. Furthermore, the argument that for large X the 
behaviour of the system can be described by an effective dichotomous process, is true for a general 
potential as well. 

Let us discuss the relevance of our results for the motion of molecular motors. In [1C] it was argued 
that the basic features of the movement can be described by a sum of AT — 1 dichotomous processes, 
where N > 3 is the number of conformational changes of the motor protein. Even if one takes a simple 
sawtooth ratchet as interaction potential, this model leads to a current reversal for a proper choice of the 
parameters. According to Astumian and Bier this means that, after redimensionalizing the equation of 
motion (2T), two proteins with a slight different geometry — and thus different friction coefficients — 
may drift in opposite directions [15]. 

On the other hand, our results emphasize that rate and direction of the movement are very sensitive to 
the details of the potential. They suggest that slight modifications of the binding energy profile which the 
motor protein feels when it walks on the periodic polymer may change its direction of movement. Even 
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though the polymer tracks are much more symmetry broken than the potentials that we have investigated, 
it should be clear that a sawtooth ratchet is a too rough approximation for the actual interaction potential. 
Due to our results, one would expect that a small deviation of the simple sawtooth ratchet has a large 
effect on the current, especially in the region of intermediate x. Therefore one should be careful: A 
quantitative agreement between calculated currents and experimental findings may be accidental. But the 
qualitative aspect of the model may be correct. 

Clearly, the model is too simple to provide a realistic description of the situation in a cell. The kinetics 
of the ATP hydrolysis in the cell is very sensitive to temperature changes of the system. For this reason, 
X and T cannot be chosen independently. Another serious difficulty is that the forces exerted by the 
protein on the track during the ATP hydrolysis change not only the barrier heights but also the shape of 
the interaction potential. This conformational flexibility is not taken into account in the model, where the 
interval lengths A, and the relative heights of the potential teeth are time independent. But a more realistic 
model will contain many free parameters and, as we have shown, the results may depend strongly on each 
of these parameters. Therefore a more detailed experimental knowledge of the system is required. 
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